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Important tasks in cavity quantum electrodynamics include the generation and control of quantum 
states of spatially-separated particles distributed in different cavities. An interesting question in 
this context is how to prepare entanglement among particles located in different cavities, which 
are important for large-scale quantum information processing. We here consider a multi-cavity 
system where cavities are coupled to a superconducting (SC) qubit and each cavity hosts many 
SC qubits. We show that all intra-cavity SC qubits plus the coupler SC qubit can be prepared 
in an entangled Greenberger-Horne-Zeilinger (GHZ) state, by using a single operation and without 
the need of measurements. The GHZ state is created without exciting the cavity modes; thus 
greatly suppressing the decoherence caused by the cavity-photon decay and the effect of unwanted 
inter-cavity crosstalk on the operation. We also introduce two simple methods for entangling the 
intra-cavity SC qubits in a GHZ state. As an example, our numerical simulations show that it is 
feasible, with current circuit-QED technology, to prepare high-fidelity GHZ states, for up to nine 
SC qubits by using SC qubits distributed in two cavities. This proposal can in principle be used 
to implement a GHZ state for an arbitrary number of SC qubits distributed in multiple cavities. 
The proposal is quite general and can be applied to a wide range of physical systems, with the 
intra-cavity qubits being either atoms, NV centers, quantum dots, or various SC qubits. 

PACS numbers: 03.67.Bg, 42.50.Dv, 85.25.Cp, 76.30.Mi 


I. INTRODUCTION 

Superconducting devices can be fabricated using modern integrated circuit technology, their properties can be 
characterized and adjusted in situ, and their coherence time has recently been significantly increased [1-9]. Moreover, 
various single- and multiple-qubit operations with state readout have been demonstrated [10-15], and nonlinear optical 
processes in a superconducting quantum circuit have been investigated [16]. In addition, Circuit QED, consisting of 
microwave resonators and SC qubits, is particularly attractive and considered as one of the leading candidates for 
QIP [17-23]. The strong and ultrastrong coupling between a microwave cavity and SC qubits has been demonstrated 
in experiments (e.g., [24-26]). In addition, using SC qubits coupled to a single cavity or resonator (hereafter, the 
terms cavity and resonator are used interchangeably), a number of theoretical proposals have been presented for 
realizing quantum gates and entanglement [17-19,27-32], and two- and three-qubit quantum gates and three-qubit 
entanglement have been experimentally demonstrated [33-37]. 

In recent years, there is much interest in large-scale QIP, which usually involves many qubits. Note that placing all 
qubits in a single cavity may cause many problems, such as increasing the cavity decay rate and decreasing the qubit- 
cavity coupling strength. Therefore, for cavity or circuit QED-based large-scale QIP, the qubits should be distributed 
in different cavities, and the ability to perform nonlocal quantum operations for these qubits is a prerequisite to realize 
distributed quantum computation. During the past few years, attention has been paid to the preparation of entangled 
states of two or more cavities, or of qubits located in different cavities, and implementation of quantum logic gates 
on photons/qubits distributed over different cavities in a network. Specifically, rapid progress has been achieved in 
the following two directions: 

(i) Manipulating and generating nonclassical microwave field states with photons distributed in different cavities. 
By using a SC qubit (artificial atom) coupled to cavities, schemes have been proposed for synthesizing different 
entangled photonic states of two SC resonators [38], and for generating multi-particle entangled states of photons in 
different cavities [39,40]. By employing the idea of Ref. [41] the so-called NOON state of photons in two resonators 
has been experimentally created [42]. In addition, how to perform quantum logic operations on photons located in 
different cavities has been investigated [43]. 

(ii) Quantum state engineering and quantum operations with qubits distributed in different cavities. By using a 
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SC qubit to couple two or more cavities/resonators, proposals have been presented for generating GHZ states with 
multiple SC qubits coupled to multiple resonators via employing cavity photons and through step-by-step control 
[40,44], and for quantum information transfer between two spatially-separated SC qubits distributed in two cavities 
[45]. Recently, GHZ states of three SC qubits in circuits consisting of two resonators have been experimentally 
prepared [7]. Furthermore, using an intermediate SC qubit coupled to two planar resonators, quantum teleportation 
between two distant SC qubits has recently been demonstrated in experiments [46]. 

GHZ states are not only of great interest for fundamental tests of quantum mechanics [47], but also have applications 
in QIP [48], quantum communications [49], error-correction protocols [50], quantum metrology [51], and high-precision 
spectroscopy [52]. During the past years, experimental realizations of GHZ states with eight photons using linear 
optical devices [53,54], fourteen ions [55], three SC qubits in circuit QED [7], five SC qubits via capacitance coupling 
[56], and three qubits in NMR [57] have been reported. Theoretically, proposals for generating entangled states with 
SC qubit circuits have been presented [58]. In addition, based on cavity QED or circuit QED, a large number of 
theoretical methods have been presented for creating multi-qubit GHZ states with various physical systems (e.g., 
atoms, quantum dots, and SC devices) that are coupled to a single cavity/resonator mode [59-68]. However, we note 
that how to generate GHZ states with qubits in different cavities has not been thoroughly investigated. 

Motivated by the above, here we present an efficient method to entangle SC qubits in a multi-cavity system, where 
cavities are coupled to a SC qubit and each cavity hosts many SC qubits. We show that the cavity-induced effective 
conditional dynamics between the intra-cavity SC qubits and the coupler SC qubit can be employed to entangle all 
the SC qubits in a GHZ state. In this work, we also introduce two simple methods for entangling the intra-cavity SC 
qubits in a GHZ state. As an example, our numerical simulations show that it is feasible, with current circuit-QED 
technology, to prepare high-fidelity GHZ states, for up to nine SC qubits by using SC qubits embedded in two cavities. 
To the best of our knowledge, based on circuit QED, the experimental demonstration of GHZ states has only been 
reported for three SC qubits [7,34]. 

This proposal has the following advantages: (i) The GHZ state preparation does not require step-by-step control, 
which involves only one operation for entangling all qubits and a few basic operations for entangling the intra-cavity 
qubits; (ii) The entanglement is prepared without exciting the cavity photons, and thus the decoherence induced by 
cavity decay and the effect of unwanted inter-cavity crosstalk on the operation are greatly suppressed; (iii) Because 
none of the intra-cavity qubits is excited during the operation, decoherence from the qubits is much reduced; (iv) 
More interestingly, this proposal can in principle be used to implement a GHZ state for an arbitrary number of 
qubits distributed in multiple cavities by using a single coupler qubit, which is important for the future realization 
of large-scale QIP; and (v) We further stress that this proposal is quite general and can be used for other kinds of 
qubits, such as atoms, NV centers, and quantum dots. 

This paper is organized as follows. In Sec. 2, we introduce the physical model considered in this work and derive 
the effective Hamiltonian used for the entanglement production. In Sec. 3, we show how to generate GHZ states for 
all the intra-cavity SC qubits and the coupler SC qubit based on the effective Hamiltonian. In Sec. 4, we further 
introduce two simple methods for generating GHZ states of intra-cavity SC qubits. In Sec. 5, as an example, we 
numerically analyze the experimental feasibility of preparing a GHZ state of up to nine SC qubits with SC qubits 
distributed in two cavities. A concluding summary is given in Sec. 6. For the numerical calculations, here we use the 
QuTiP software [69,70]. 


II. PHYSICAL MODEL AND EFFECTIVE HAMILTONIAN 

We consider a system composed of N cavities and assume that cavity j hosts m,j SC qubits denoted as ji, j 2 ,..., 
and jm r These cavities (1,2,..., N) are coupled to a common SC qubit A (coupler qubit), as shown in Fig. 1. Each 
qubit considered here has three levels, which are denoted as | g ), | g '), and |e) (Fig. 2). The two logical states of each 
intra-cavity qubit are represented by the two levels | g) and | g') 1 while those of the coupler qubit are represented by the 
two levels | g') and |e). The third level |e) for each intra-cavity qubit or |g) for the coupler qubit acts as an auxiliary 
level for realizing, e.g., a conditional phase shift. The level spacing of the coupler qubit A is different from those of 
the intra-cavity qubits. The | g) o |e) transition of the intra-cavity qubits is coupled to their respective cavities with 
coupling strength /i, while the |<?) -SA |e) transition of the coupler qubit A is coupled to all the cavities with coupling 
strength /ia (Fig- 2). We here assume that the level | g') of each qubit is not affected during the operation, which 
can be a good approximation when the transitions between | g') and other levels are sufficiently weak or the relevant 
transition frequencies are highly detuned from the cavity frequency. The level spacings of superconducting qubits can 
be rapidly adjusted by varying the external control parameters (e.g., the magnetic flux applied to a superconducting 
loop for phase, transmon or flux qubits; see e.g. [20,71-73]), so that the transitions associated with | g') can be tuned 
far off-resonance with the resonators. In the interaction picture with respect to the free Hamiltonian of the system 
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FIG. 1: (Color online) Diagram of a coupler qubit A (circle at the center) and N cavities (1, 2, N) each hosting qubits. Each 
cavity here is a one-dimensional coplanar waveguide transmission line resonator. The circle A represents a SC qubit, which is 
capacitively coupled to each cavity. A dark dot represents an intra-cavity SC qubit. For simplicity, only three qubits are drawn 
in each cavity. 



(a) 


(b) 


FIG. 2: (Color online) (a) Each cavity is dispersively coupled to the |g) |e) transition of the coupler qubit A, with coupling 

strength g A and detuning A A = u> c — co eg ,A- Here, us c is the cavity frequency while tu eg ,A is the | g) \e) transition frequency 
of the coupler qubit A. (b) Each cavity is dispersively coupled to its hosting intra-cavity qubits’ | g) |e) transition, with 
coupling strength g and detuning A = oj c — uj eg . Here, oj eg is the | g) |e) transition frequency of the intra-cavity qubits. In 
addition, <5 = A — A a- In (a) and (b), a horizontal solid line represents an energy level of a qubit; and each vertical (red) 
double-arrow line, linked to a middle solid line and a top (red) dashed line, represents the cavity frequency uj c . For simplicity, 
we here consider the case when the level spacing between the two lowest levels is smaller than that between the upper two 
levels. This type of level structure is available in superconducting charge qubits or flux qubits. Alternatively, the level spacing 
between the two lowest levels can be larger than that between the upper two levels, which applies to superconducting phase, 
transmon or Xmon qubits. In (a) and (b), the ground level is labeled by \g' A ) (|</)) and the first excited level is denoted as \gA ) 
(|g)). One could also denote the ground level as \g A ) (|s)) and the first excited level as | g' A ) (!</))■ 


(not shown for simplicity), the Hamiltonian is given by 

N N 

H = + e iAt a]jj^j + ^ g A (eT lAAt a 3 a\ + e lAAt a]a A ^j , (1) 

i=i j =i 

where jf = <J+, Jj = cr 3z , a+ = \e) u (g \, a 3z = \g) h (e \, a\ = \e) A {g\, a A = \g) A (e|, A = w c - w eg , and 

Aa = oj c — w e g,A- Here, uj c is the cavity frequency; and cu eg and u> e g,A are the | g) -o- \e) transition frequencies for the 
intra-cavity qubits and the coupler qubit A , respectively (Fig. 2). 

Under the large-detuning condition A, A a , g, g A , |<5| = | A — Aa| , the dynamics governed by H is equivalent to 
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that decided by the following effective Hamiltonian [74-76] 

i/eff = -y (Gj a + a j — Ejajaj ) 

i=i 
2 JV 

+X7 J2 ( Ga s + ) 

A 3 =1 

2 W 
3 =i 
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-X'£(e~ iSt J j a++H.c.) 

3 =1 
N 

+ Ajfc (aja'l + H.c .) (Ea — Ga) j (2) 
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where S = A- A a , Gj = J2?Ji \g) jt id I» E j = El =*i l e )j 4 ( e l. Ga = |s)a (ffl > e a = \e) A (e|, A jfc = [i\/A a , and 
A = ^ . The terms in the first and second lines of (2) account for the ac-Stark shifts of the level |g) (|e)) of 

the intra-cavity qubits and the coupler qubit A induced by the corresponding cavity modes, respectively. In addition, 
the terms in the third line describe the effective dipole-dipole interaction between the intra-cavity qubits located in 
the same cavities. The terms in the fourth line describe the dipole-dipole interaction between the intra-cavity qubits 
and the coupler qubit A, and the terms in the last (fifth) line characterize the coupling between any two cavities. 

If each cavity is initially in a vacuum state, the Hamiltonian (2) reduces to 


= -iE £ i- 
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( 3 ) 


According to [74], if |<5| » A, the terms in the last line can be replaced by ^-Ej=i ^j°A>X^i J^cta]- 

Thus, the dynamics governed by Hamiltonian (3) is approximately equivalent to that by the following Hamiltonian 
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which can be rewritten as 
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(5) 


When the level |g) of the coupler qubit and the level |e) of the intra-cavity qubits are not populated, the Hamiltonian 
(5) reduces to 


A, = -M EA + ^E A -£a, 


N 


A a"" ' 6 


Nii 2 \ 2 N mj 

I e ) A ( e l + -J I e) A (e| ® EE l5>j 4 (si • 
A j =1 *=1 


( 6 ) 


This effective Hamiltonian can be turned off by tuning the qubit levels in such a way that the transitions between these 
levels are highly off-resonant with the cavities, and hence the qubits are effectively decoupled from the corresponding 
cavities. 


III. ENTANGLING INTRA-CAVITY QUBITS AND THE COUPLER QUBIT 


Let us go back to the setup in Fig. 1. Initially, the qubit system is decoupled from the cavity system, and each 
cavity is in the vacuum state. Assume that each intra-cavity SC qubit is in the state |+) = l/-\/2(|</) + | g)) and the 
coupler SC qubit is in the state a \ g') A + /3 \e) A (|a| 2 + |/3| 2 = 1). These states can be prepared from the qubit ground 
state with classical pulses. For simplicity, let us consider the case of | g) being the ground state (a case applied to the 
flux qubits considered in Sec. 5, with three levels illustrated in Fig. 6). The state |<?) can be transformed to |+) by 


applying a 7r/2 pulse tuned to the | g) A |gj transition [29]. The preparation of the state a\g') A + (3 \e) A consists 

of two steps: (i) apply a n pulse, tuned to the \g) A -A \e) A transition, to flip the state | g) A to le)^; (ii) employ a 
classical pulse to drive the | g') A A le)^ transition, with the Rabi frequency fl and duration t satisfying a = cos (fit.) 
and /3 = sin (fit.). 

The initial state of the whole qubit system is thus given by 


N rrij 

(a | g') A + /? \e) A ) <g> \+) k ■ 

j=i i=i 


(7) 


Now adjust the level spacings of qubits, so that the qubit-resonator coupling is turned on and the dynamics of the 
qubit system is governed by the effective Hamiltonian (6). One can see that under the Hamiltonian (6), the state (7) 
evolves into 


N 271 j N 271 j 
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(8) 
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FIG. 3: (Color online) Sequence of operations on each qubit (from left to right). Here, t satisfies Eqs.(9) and (10), which is 
the qubit-cavity interaction time required for producing the GHZ state; while td (within 1-3 ns [77,78]) is the time required 
to adjust the qubit level spacings. Note that the level spacings of qubits are tuned simultaneously. For simplicity, here and in 
Figs. (4) and (5), we assume that the time needed for adjusting the level spacings is the same for both intra-cavity qubits and 
coupler qubit. 


where we have used H' eff \g') = 0. Setting 


M, 

Aa 

A 2 

iV 


where m is an integer, Eq. (8) can be expressed as 


N mj 

Jim 

3= 1 i= 1 


2m7r, 

(9) 

7T, 

(10) 

N mj 

,+^nni-u 

(ID 


3= 1 i=l 


where |-) = l/\/2(| g') - I#)). Since |-) is orthogonal to |+), the state (11) is a multi-particle GHZ entangled state 
for the coupler qubit and the qubits distributed in multiple cavities. One can see that the entangled state preparation 
here is based on a n -phase shift on the state | g) of each intracavity qubit conditional upon the coupler qubit being 
in the state \e) A . Note that by applying a classical pulse to the coupler qubit, the states | g') A and \e) A can be easily 
converted into the states 1+)^ and |—) A , respectively. The operation sequence for preparing the GHZ state (11) is 
illustrated in Fig. 3. 

Note that when the coupling of qubits to the mode in a cavity is spatially dependent, then different intra-cavity 
qubits will acquire different conditional phases. The coupler qubit would acquire a single-qubit phase, when its 
coupling to the cavity mode is deviated from the preset value. To eliminate the effect of this phase, the interaction 
time should be adjusted so that this phase is equal to 2mn 1 with m being an integer. 

It should be mentioned that, in order to maintain the initial states and the prepared GHZ states of the qubit 
system, the coupler qubit and the intra-cavity qubits should be decoupled from their respective cavities before/after 
the entanglement production, which requires the qubit-cavity coupling to be switchable. This requirement can be 
readily achieved, by prior adjustment of the level spacings of the qubits [21,71-73] or the frequencies of the cavities. 
We note that the rapid tuning of microwave cavity frequencies has been experimentally demonstrated (e.g., in less 
than a few nanoseconds for a superconducting transmission line resonator [77,78]). 

By preparing the initial state of the coupler qubit A with different values of a and /3, the degree of entanglement 
for the GHZ state (11) can be adjusted and thus this protocol can be used to generate GHZ entangled states with 
an arbitrary degree of entanglement. As shown above, during the entanglement preparation no photons are excited in 
each cavity, no measurement is needed, and only a single-step operation is required. 


IV. ENTANGLING INTRA-CAVITY QUBITS 

In this section, we will briefly introduce two methods for preparing all intra-cavity SC qubits in a GHZ state. The 
first method requires a measurement on the coupler SC qubit, while the second one does not need any measurement. 


A. Method 1 
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FIG. 4: (Color online) (a) Sequence of operations on each intra-cavity qubit, (b) Sequence of operations on the coupler qubit. 
The green curves in (a) and (b) correspond to the operation for producing the GHZ state (11). In (b), the 7t/2 pulse transforms 
the state \g') A —» (| g') A + |e) A ) /\/2 and the state |e ) A (| g') A — |e) A ) /y/2, and the square box with a meter represents a 
measurement on the coupler qubit along the basis {| g') A , |e)^}. Note that the operations in (a) and (b) are performed from 
left to right. 


All qubits including the coupler qubit are first prepared in the GHZ state (11). One can see from Eq. (11) 
that through a unitary tranformation | g') A i\g') A + \e) A ) /\/2 and \e) A —> (\g') A — |e) A ) /\/2 and then a 
measurement on the coupler qubit A, the intra-cavity qubits will be prepared in one of the two GHZ states 
|GHZ) ± = nl=i l+)ji /3nyLi n(=i | — ) J - i (depending on the measurement outcome). The operation se¬ 

quence for preparing the states IGHZ^ is illustrated in Fig. 4. 

The two GHZ states |GHZ) + and |GHZ) _ can be converted into each other through the local operation on any one 
of intra-cavity qubits (say qubit j i): |—) ^ —> — |—) ^ and |+) ^ —> |+) J1 - I 11 this sense, the intra-cavity qubits can 
be prepared in a GHZ entangled state deterministically. 

As discussed here, a measurement on the coupler qubit is necessary to prepare the intra-cavity qubits in a GHZ 
state. Note that fast and highly accurate measurements on the state of a SC qubit are experimentally available at this 
time (e.g., see [14]). In the following, we will propose an alternative approach for entangling the intra-cavity qubits, 
which does not require any measurement. 


B. Method 2 

Assume that one intra-cavity SC qubit, say qubit li in cavity 1, is initially in the state |(/) , each of all remaining 
intra-cavity SC qubits is initially in the state |+), and the coupler SC qubit A is in the state a | g') A + /3 |e) A . Then 
the initial state of the whole system is thus given by 

mi N 171 j 

I initial) = II l+)i, ® IIII l+)j, ® ( a 15% + Z 3 \ e ) a ) ■ ( 12 ) 

i= 2 j= 2 i= 1 

The procedure for preparing intra-cavity qubits in a GHZ entangled state is listed as follows: 

Step 1: Keep qubit li decoupled from cavity 1, while adjust the level spacings of other qubits such that their 
dynamics is governed by the Hamiltonian (6) (not including qubit li) for an interaction time t satisfying Eqs. (9) and 
(10). By a similar derivation as shown in Eq. (8), one can easily find that the state (12) changes to 

( mi N rrij mi N rrij 

«n\ + h ® n n i+>* ® i a') a + pH h>i, ® n n i->* ® \^ A 

t=2 i=l i=2 j=2i=l 

(13) 



Then, adjust the level spacings of the qubits such that the qubit system is decoupled from the cavities. 

Step 2: Perform the operations, \g')i 1 —> |g)j and | g') A —> \g) A , by applying classical pulses to qubit li and 
the coupler qubit. In addition, perform a swap operation |< 7 ) li \e) A — > —i |e) 1 | g) A [75], which can be achieved by 

















Step 1 



Step 2 — 

tswap 



FIG. 5: (Color online) (a) Sequence of operations on qubit li. (b) Sequence of operations on the intra-cavity qubits (except 
for qubit li). (c) Sequence of operations on the coupler qubit A. The green curves in (b) and (c) correspond to the operation 
described by step 1, which is for preparing the state (13). The red curves in (a) and (c) correspond to the operation of step 2 
for preparing the state (14). The n pulse in (a) transforms the state |</ , ) 1 —» |<?) ll , while the n pulse in (c) transforms the state 
| g') A —> | g) A . The right-hand square curves in (a) and (c) represent a swap operation described by |g) 1 \e) A —» —i \e) 1 \g) A 
[75]. In (a) and (c), t SW a P is the swap operation time. Note that the operations in (a), (b), and (c) are performed from left to 
right. 


adjusting the level spacings of qubit li and the coupler qubit, such that the transitions | g) —> |e) of qubits li and the 
coupler qubit are dispersively coupled to cavity 1 with the same detuning. Then the state (13) becomes 

( mi N rrij rm N mj \ 

«i5)i! n i+>i 4 ® n n i+>* - ^ n i->l ® n n i->* > 

i =2 i=2 *=1 i =2 j—2 i—1 J 

(14) 

which shows that all intra-cavity qubits are deterministically prepared in a GHZ entangled state and disentangled 
from the coupler qubit. Since no measurement is involved, the operations above are unitary. The operation sequence 
for preparing the GHZ state is illustrated in Fig. 5. 

From the description given above, one can see that the entanglement production does not employ cavity photons. 
In addition, the GHZ state preparation here does not depend on the number of intra-cavity qubits, which requires 
only a few basic operations. Hence, the methods presented here for entangling the intra-cavity qubits are quite simple. 


V. EXPERIMENTAL FEASIBILITY OF ENTANGLING MULTIPLE QUBITS: AN EXAMPLE 

To illustrate the experimental feasibility of our scheme, we consider a system of two cavites (i.e., two one-dimensional 
transmission line resonators) each hosting M superconducting flux qubits (with M ^ 4) and coupled by a supercon¬ 
ducting flux qubit A. Figure 6 shows the setup for each cavity hosting four flux qubits. The three levels |g), | g '), 
and |e) of each flux qubit are depicted in Fig. 7. The | g) -n- \g') transition of each flux qubit can be made weak 
by increasing the potential barrier and thus its coupling with the cavities is negligible. In reality, the | g') -n- |e) 
transition needs to be considered because the coupling between this transition and each cavity may turn out to affect 
the operation fidelity. Note that the | g') -n- |e) transition is much weaker than the | g) -o- |e) transition due to the 
potential barrier between the two wells. 

When the unwanted coupling of the | g') o |e) transition with the cavities and the unwanted inter-cavity crosstalk 
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FIG. 6: (Color online) (a) Setup of two one-dimensional transmission line resonators each hosting four flux qubits (dark dots) 
and coupled to flux qubit A (in the middle circle). The coupler qubit A is connected to the two resonators via capacitors C\ and 
C 2 , respectively. Each qubit could be a radio-frequency superconducting quantum interference device (rf SQUID) consisting of 
one Josephson junction enclosed by a superconducting loop as depicted in (b), or a superconducting device with three Josephson 
junctions enclosed by a superconducting loop as shown in (c). Ej is the Josephson junction energy (0.6 < a < 0.8). The 
superconducting loop of each qubit, which is a large circle for (b) while a large square for (c), is located in the plane of the 
resonators between the two lateral ground planes. Each intra-cavity qubit is coupled to its cavity via the magnetic flux through 
the superconducting loop of each qubit, which is created by the cavity magnetic field threading the superconducting loop. 
The intra-cavity qubits are placed at locations where the cavity magnetic fields are the same to achieve an identical coupling 
strength for each qubit. For each qubit, a tunable-coupler dc current line, e.g., the red dashed line in (b) or (c) placed on the 
qubit loop, creates a dc magnetic pulse threading the loop of each qubit, which is used for tuning the qubit level spacings. Note 
that the qubit level spacings can be tuned by varying the magnitude of the dc magnetic pulse through changing the current 
intensity. In addition, for each qubit, a microwave bias ac current line, e.g., the dark dashed line in (b) or (c) on the qubit loop, 
creates an ac magnetic pulse threading the loop of each qubit, which is used to prepare the initial state of each qubit or/and 
manipulate the state of each qubit during the GHZ state preparation. 




FIG. 7: (Color online) Illustration of the qubit-cavity dispersive interaction. The left is for the coupler flux qubit A , while the 
right is for the intra-cavity flux qubits. The tunneling between the two lowest levels is made weak by increasing the potential 
barrier, such that the level | g') can be stored for a long time. 
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are included, the Hamiltonian (1) is modified as 

2 2 

H = ^2n(e- iAt aj Jf+R.c.)+^2n A ( e - iAAt ci j a++R.c.) 

3 =1 3 =1 

2 2 
+ p (e~ lAt aj Jd + H.c.^ + ^ p A ^ e _ * A ' 4t a.,<4 + H.c.^ 
i=i i=i 

+/ii 2 (aiaj +H.c.), (15) 

where J+ = Ej=i E^i <4> 4 + = E./=i Eti = le)^ (p'|, and d+ = |e) A (g'l. The terms in the first (second) 

pair of parentheses in the second line describe the unwanted coupling between the | g') |e) transition of intra¬ 

cavity qubits (coupler qubit) and the respective cavities with coupling strength p ( p A ) and detuning A = w c — LO eg > 
(A a = lj c — ut eg ' >A )• The terms in the last pair of parentheses describe the inter-cavity crosstalk between two cavities, 
with the inter-cavity coupling constant pvi- 

The dynamics of the lossy system, with finite qubit relaxation, dephasing and photon lifetime being included, is 
determined by the following master equation 

2 

+ ^ ( K jP \ a i\ T P + Qi (16) 

j =i 

with 

2 M 

p = J2J2{^ g P-[ a ji}+'yeg'C[a ji ]+'y gl gC[a ji }} 

3= 1 i=l 

2 M 

+ {le,tp {PjipPji ~ PjiPl^ 1 ~ PPji/ty} 

j=1 i=l 

2 M 

+ E E Kb* (<4^4 - G'.p/2 - pG'jJ 2)} , (17) 

i=i i=i 

<3 = 7e 9 ,A/3 [cr“] +7eg/,2l£ [?l] + 7 9 's,a£ [ct a ] 

+7e, v ,A ( E A pE A - E A p/2 - pE A /2) 

+7 g',ip, A ( G' A pG' A — G' A p/2 - pG' A /2 ), (18) 

where ct ,- 4 = |p 4 (ff'| ,aA = |s) A (. g '\; = \e) h (e|, G'. = \g') k {g '\, G A = \g') A (g '\; and £ [A] = ApA+ - A+Ap/2 - 

pA + A/2, with A = aj, aj i ,aj i ,aj i , a A ,a A ,cr A . Here, Kj is the photon decay rate of cavity a g . In addition, 7 g > g is the 
energy relaxation rate of the level \g') of qubits, 7 eg (' y eg >) is the energy relaxation rate of the level \e) of intra-cavity 
qubits for the decay path |e) —> |g) (|</)), and 7 e>v {l g ',g>) is the clephasing rate of the level |e) (|</)) of intra-cavity 
qubits. The symbols 7 g i gtA , 'Yeg,A, leg',A, 7 e,g>,A, and 7 g ', Vt A denote the corresponding decoherence rates of the coupler 
qubit A. 

The fidelity of the operation is given by [79] 


ap 

It. 


H, 


P = 4(V’id|p|V’id), (19) 

where \ipid) is the output state of an ideal system (i.e., without dissipation, deplrasing, and crosstalks) as discussed 
in the previous section, and p is the final density operator of the system when the operation is performed in a 
realistic physical system. Without loss of generality, consider now that the qubit system is initially in the state 
n] = i Tlf=i 1+4 (| g') a + \ e ) a) /\/2 (with M < 4), and the two cavities are initially in a vacuum state, for which the 
ideal state |?/>id) is the one given in Eq. (11) with a = /? = l/\/2. 

We consider ident ical intra-cavity flux qubits. Given A, A^, and to, the coupling constant p is determined by 
p = \/2N |<5| A A /mA/ (A + A^) [derived from Eqs. (9) and (10)]. We set p = 0.1 p and p A = 0.1 p A , which is a 
good approximation by increasing the potential barrier such that the transition matrix element between the two levels 
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FIG. 8: (Color online) Fidelity versus ha/h- Here, n is the total number of qubits including the coupler qubit. Also n = 3, 5, 7 
and 9 correspond to the cases of each cavity hosting one, two, three and four qubits, respectively. To have an achievable good 
fidelity with reasonable parameters, we increase A/(27r) = 2 from 1 to 2 GHz, when going from 5 to 7 qubits. For n = 3, 5, 
7 and 9, a high fidelity ~ 97.2%, 94.3%, 91.3%, and 90.1% can be achieved with /za/m being 1.516, 1.641, 1.174, and 1.242, 
respectively. 


| g') and |e) is smaller than that between the two levels | g) and |e) by one order of magnitude (Fig. 7). We choose 
A/ (27t) = 1.0 GHz for n = 3 and n = 5, while A/ (27r) = 2.0 GHz for n = 7 and n = 9. Here, n = 2 M + 1 is the total 
number of qubits including the coupler qubit. In addition, we set A = A + 27T x 1.5 GHz and A a = Aa + Szt x 1.5 GHz. 
We now choose nj 1 = 15 /zs, 7 “* = 5 /zs, 7 ^ = 7.5 /us, 7 ~g = 5 /us, 7 ~ g , = 7.5 /us, and 7 ^ = 15 /zs (a conservative 

consideration, e.g., see Ref. [9]). Our numerical calculations show that when /Z 12 is smaller than /z a by two orders 
of magnitude, the effect of the inter-cavity crosstalk on the operation is negligible. Thus, we now set /U 12 = 0.01/za- 
With the parameters chosen here and by numerically optimizing the parameters m, /u, and Aq, the fidelity versus 
/za//z is plotted in Fig. 8 for n = 3, 5, 7 and 9. From Fig. 8 , one can see that for n = 3, 5, 7 and 9, a high fidelity 
~ 97.2%, 94.3%, 91.3%, and 90.1% can be achieved with the optimized values of ha/I u being 1.516, 1.641, 1.174, and 
1.242, respectively. We remark that the fidelity can be further increased by improving the system parameters. 

Figure 8 shows that with the parameter values chosen above the large detuning condition is not well satisfied for the 
optimized fidelity, e.g., Aa/Ha ~ 3.95 for n = 9, which implies that the state evolution determined by the effective 
Hamiltonian ( 6 ) can be a good approximation with a suitable choice of parameters, even when the qubit-cavity 
detunings are not much larger than the coupling strengths. This result has a quantitative explanation. Beyond the 
large detuning regime, when the coupler qubit is initially in the state |e), the total cavity-qubit system undergoes 
Rabi oscillations in the corresponding single-excitation subspace. The associated Rabi frequencies have a dependence 
on the number of the intracavity qubits in the state | g). With a suitable choice of the qubit-cavity coupling strengths 
and detunings, all of the state components with the coupler qubit being initially in |e) can return to their initial forms 
almost at the same time, with the resulting phase shift being related to the corresponding Rabi frequency. 

As discussed in [39,80], the condition H 12 < 0.01 ha can be met with the typical capacitive cavity-qubit coupling. 
Figure 8 shows that at the optimum points, the coupling strengths are {/z/ 27 t, /za/ 27 t} ~ {180.4 MHz, 273.5 MHz} 
(n = 3), {156.4 MHz, 256.75 MHz} (n = 5), {232.7 MHz, 273.2 MHz} {n = 7), and {269.0 MHz, 334.0 MHz} (n = 9). 
The coupling strengths of these values are readily achievable in experiments because a coupling strength ~ 636 MHz 
has been reported for a superconducting flux device coupled to a one-dimensional transmission line resonator [26]. 
For a flux qubit, the | g) 0 |e) transition frequency could be between 5 GHz and 20 GHz. Thus, we can choose 
w c /27t ~ 7.5 GHz. For the value of nj 1 used in the numerical calculation, the required quality factor for each cavity 
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TABLE I: Qubits li, I2, I3, and I4 are four qubits placed in one cavity; while Qubits 2i, 22 , 23 , and 24 are four qubits placed 
in the other cavity. Here p takes the same value used in Fig. 8 . 


TABLE II: Here A takes the same value used in Fig. 8 . 


is Q = KjUJ c ~ 7.1 x 10 5 , which is available in experiments according to previous reports [81,82]. Therefore, the 
high-fidelity creation of GHZ states of up to nine qubits by using this proposal is feasible with current circuit QED 
technology. 

To see how well this method works in a more realistic situation, we now consider inhomogeneous coupling of qubits 
to the mode in each cavity, non-uniform distribution of qubit frequencies, imperfect preparation of the initial states, 
error in the operation time, and the existence of thermal photons in each cavity. The density operator of a thermal 


state of each cavity is described by p 


E 


n=0 


n n 

(l+n)" + 1 


| n) ( n \, with n being the average photon number and |n) being 


an n-photon state. In our numerical simulation, we choose n = 0.1. The initial state of the qubit system is modified 
as n-=i nf=i (yj ^ 2 ^ W) ji + \J I <7)jq) (\J 1 9') a + \J |e) A ) • F° r simplicity, we here consider the identical 
error e for the preparation of the initial state of each qubit. Without loss of generality, we will numerically investigate 
how the maximum fidelity in each subfigure of Fig. 8 is affected by the above-mentioned factors. 

As mentioned above, the coupling constant pa of the coupler qubit, corresponding to the maximum fidelity in each 
subfigure of Fig. 8 , is calculated to be Pa = 1.516p for n = 3; pa = 1-641 p for n = 5; pa = 1.174 p for n = 7; and 
PA = 1.242/i for n = 9. Here, the value of p is shown in Fig. 8 . In addition, the detuning A a for the coupler qubit 
takes the same value as in Fig. 8 . 

The coupling constants and the detunings for the intra-cavity qubits are listed in Tables 1 and 2, where up to 1% 
inhomogeneous coupling constants and up to 5% non-uniform detunings are considered. In Tables 1 and 2, the values 
of p and A are the same as those in Fig. 8 . 

With the parameters chosen above, in Fig. 9 we present a numerical simulation of the fidelity versus t/T for e = 0, 
0.05, and 0.1. Here, t is the operation time, while T is the optimal operation time corresponding to the maximum 
fidelities in Fig. 8 , which are 46.55 ns, 67.94 ns, 167.19 ns, 106.43 ns for n = 3,5,7, and 9, respectively. From Fig. 
9, one can see that the fidelity is insensitive to the error e but is significantly affected by the error in operation time. 
For t = 0.95T or 1.05T (i.e., 5% operational time error), the fidelity drops down to < 0.5. Note that for t = T, good 
fidelities > 0.939 for n = 3, > 0.845 for n = 5, > 0.778 for n = 7, and > 0.760 for n = 9 can be obtained. 

It is worthwhile to discuss the advantage of utilizing negative detunings versus positive detunings. For the flux 
qubits with three levels | g), | g'), and |e) shown in Fig. 7, the purpose of using a negative detuning is to increase 
the detuning of the |g') O |e) transition frequency from the cavity frequency, in order to reduce the effect of this 
unwanted transition on the operation fidelity. It can be seen from Fig. 7 that the detuning of the | g') o |e) transition 
frequency from the cavity frequency would be smaller when using a positive detuning (i.e, the case when the cavity 



13 



FIG. 9: (Color online) Fidelity versus t/T. Here, t is the operation time, while T is the operation time corresponding to the 
maximum fidelities in Fig. 5, which are 46.55 ns, 67.94 ns, 167.19 ns, 106.43 ns for n = 3,5,7, and 9, respectively. The blue, 
green, and red curves correspond to e = 0,0.05, and 0.1. For t = T, the fidelities corresponding to e = {0,0.05,0.1} are: 
{0.942, 0.941, 0.939} for n = 3; {0.850, 0.848, 0.845} for n = 5; {0.790, 0.789, 0.778} for n = 7; and {0.786, 0.769, 0.760} for 
n = 9. 


frequency is smaller than the | g') O |e) transition frequency), compared to using the negative detuning. 


VI. CONCLUSION 

We have proposed a general and efficient way to entangle SC qubits in a multi-cavity system. In principle, GHZ 
states of an arbitrary number of intra-caviy qubits plus the coupler qubit can be created through a single operation 
and without any measurements. Since only virtual photon processes take place, the decoherence caused by cavity 
decay and the effect of unwanted inter-cavity crosstalk are greatly suppressed. Also, the higher energy level is not 
occupied for any intra-cavity qubit; thus decoherence from the qubits is much reduced. In addition, we have introduced 
two simple methods for entangling the intra-cavity qubits in a GHZ state. Our numerical simulations show that it 
is feasible to generate high-fidelity GHZ entangled states with up to nine SC qubits in a circuit consisting of two 
resonators. We hope this will stimulate future experimental activities. The method presented here is quite general 
and can be applied to various other physical systems. We believe that the cluster-style architecture shown in Fig. 1 
has applications in fault-tolerant code for scalable quantum computing. Multiple physical qubits in each cavity can 
be used to construct a logic qubit, as required by error-correction protocols. The system can also be used to simulate 
the dynamics of the star-type coupled spin system, where many spins are coupled to a common spin. 
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